Thermal Fluctuations and Rubber Elasticity 
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The effects of thermal elastic fluctuations in rubber materials are examined. It is shown that, 
due to an interplay with the incompressibility constraint, these fluctuations qualitatively modify the 
large-deformation stress-strain relation, compared to that of classical rubber elasticity. To leading 
order, this mechanism provides a simple and generic explanation for the peak structure of Mooney- 
Rivlin stress-strain relation, and shows a good agreement with experiments. It also leads to the 
prediction of a phonon correlation function that depends on the external deformation. 
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The term rubber (elastomer) refers to amorphous, 
essentially incompressible, solids that consist of a 
crosslinked polymer network and that can sustain large, 
reversible, shear deformations. It has long been under- 
stood that the elasticity of rubber is predominantly en- 
tropic, being associated with the suppression of the en- 
tropy of the polymer network by the imposed deforma- 
tion. As a result rubber elasticity is characterized by a 
shear modulus that is proportional to temperature. 

The classical theory of rubber elasticity P, Q , devel- 
oped by Kuhn, Wall, Flory, Treloar, and many others, 
around the 1940's, qualitatively accounts for the entropic 
nature of rubber elasticity. It is based on the crucial as- 
sumption that the junctions of polymer networks do not 
fluctuate in space, but nevertheless deform affinely with 
an imposed uniform shear strain. The entropy of the 
entire rubber network is then given by the sum of the 
entropies of each polymer chain. For a uniform shear, 
i.e., for a homogeneous, volume-preserving deformation 
A, the elastic free-energy density / is given by 



,/o = - 



V 



1 



5S = -^oTrA^A, 



(1) 



where T is the temperature, 5S is the total entropy 
change due to A, V is the volume, and /^o ~ ^BT/Sf" 
is the entropic shear modulus in d dimensions, with ^ 
being the typical mesh size of the polymer network. For 
a uniaxial shear deformation along z. 



and the classical theory predicts 



/o(A) = + \ 



(2) 



(3) 



It has long been known that the classical theory does 
not work well for large deformations 0. Its failure be- 
comes most salient in the so-called Mooney-Rivlin plot 
of the stress-strain relation, in which {df / dX) / {X — A^^) 
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FIG. 1: The Mooney-Rivlin plot of the Mooney stress 
(df /d\)/{\ — A~^) for a unia^xial shear deformation, Eq. |(2J, 
versus 1/A. While the classical theory predicts a horizontal 
line, real rubbery materials universally exhibit a pronounced 
peak feature. Squares: data from Xu and Mark 3] (unit: 10^ 
Pa); Top curve: our theory fit that incorporates lowest-order 
phonon fluctuations on scales beyond crosslink spacing, using 
Eq. 1171 . Fitting parameters: B — oo, /iq ~ 8.92, /ii = 7.10. 
Also plotted are the curves for flnite bulk moduli B, which 
are calculated using Eq. lESl and same fitting parameters 
and /ii. 

is plotted versus 1/A. While from Eq.® it is clear that 
the classical theory predicts a horizontal line at fj-o in 
such a plot, almost all rubbery materials, natural or syn- 
thetic, exhibit universal and nontrivial features (e.g., a 
peak around a compressed state) , as illustrated in Fig^ 

A number of mechanisms put forward to explain these 
features, including polymer entanglement ^[^iH, non- 
Gaussian chain statistics, irreversible effects, internal en- 
ergy effects, as well as nematic order of various types, and 
crystallization However, to date, there is no broad 

consensus on the nature of the dominant mechanism re- 
sponsible for the deviation of rubber elasticity from the 
classical theory. 

Our main result, Eq. (|17|l . plotted in Fig^ is a generic 
and simple explanation for the Mooney-Rivlin data in 
terms of incompressible phonon fluctuations of the rubber 
network. It is based on the key observation that the clas- 



2 



sical theory is an effective mean-field theory that misses 
a large entropic contribution to the free energy, due to 
thermal phonon fluctuations of the polymer network hav- 
ing wavelengths longer than ^. These corrections are in 
fact comparable to the classical elastic free-energy den- 
sity, Eq. whose scale is given by the shear modulus 
A*o ~ ksT/Sf'. To see this, we note that, by the equipar- 
tition theorem, each phonon mode contributes ksT to 
the total free energy, and that the phonon density of the 
state is given by set by the mesh size ^, which is 
the short-scale cutoff for these long-wavelength fluctua- 
tions 

It is noteworthy that these elastic fluctuations were 
previously studied by James and Guth in a phantom 
network, where interactions between neighboring poly- 
mer chains are ignored. They concluded that these fluc- 
tuations are independent of the imposed strain deforma- 
tion and therefore have no effect on the elasticity. 

However, real rubber is nearly incompressible. So, in 
addition to the global volume preserving constraint on 
the macroscopic distortion, i.e. detA = 1, the elastic 
fluctuations are constrained to preserve the local density. 
As we shall show, the interplay between this local in- 
compressibility constraint and an imposed macroscopic 
shear deformation alters the spectrum of phonon fluctu- 
ations. This leads to an entropic correction to the free 
energy that depends strongly on the imposed distortion, 
and thereby qualitatively modifies the stress-strain rela- 
tion beyond that of the classical theory. 

To go beyond classical rubber elasticity and incorpo- 
rate thermal fluctuations on scale beyond ^, we consider 
the Lagrangian |9j| of an incompressible, homogeneous, 
isotropic elastic manifold (appropriate for an amorphous 
solid, when local heterogeneities are ignored) with the 
elastic energy given by the classical theory, Eq. 



L= d^X 



ipr-(X)2-i^o {Wr{X)f 



(4) 



where r{X) is the elastic configuration of a d-dimensional 
manifold satisfying the local incompressibility constraint 



detVf{X) = 1, 



(5) 



with X the equilibrium position of the mass point in the 
undeformed reference state, and p the mass density. Now 
consider a uniformly shear-strained state characterized 
by a deformation gradient A that preserves the volume, 
i.e., detA = 1. The mass-point X, which was fluctuating 
about the position X in the unstrained state, now fluc- 
tuates around the new equilibrium position A • X in the 
strained state. 

Our task is therefore to sum over all long wave-length 
elastic fluctuations in the strained state that satisfy local 
incompressibility. However, the highly nonlinear nature 
of the constraint, Eq. |SJ), is the major impediment to the 
analysis. The key to our progress in the solution of this 



problem is the resolution of the constraint ((SJ through 
the following parameterization of r{X): 



r{X) = A • e"(-^)-^X. 



(6) 



It can be proved that f{X) is volume-preserving if and 
only if the fleld v{X) is divergenceless, i.e.. 



V • v{X) = 0, or q- 0, 



(7) 



where the second form is the equivalent transversality 
(to g) condition on v^, the Fourier transform of v{X). 
It is customary to call such a v{X) field an incompress- 
ible flow The resolution of the nonlinear constraint 
Eq. @ by the linear relation Eq. Q in terms of v paves 
the way for a systematic treatment of thermal fiuctua- 
tions in an incompressible elastomer. 

Setting V — in Eq. ©, we find = A ■ X, i.e. 
the uniformly strained reference state. For small v, the 
exponential in Eq. ^ is 



f{X) = A • + v{X) + • • • ) 



(8) 



identifying A • v{X) with a phonon field displacement 
from the uniformly-strained reference state fg. 

The canonical momentum field H{X) conjugate to 
v{X) is calculated in the standard way 0- To lowest 
order in v, we find 



n(x) = 



SL 

6v{x) 



PS ■ viX), 



(9) 



where g = A"^A is the metric tensor. 

The Hamiltonian H is related to the La- 
grangian via the Legendre transformation: 
H[ii,v] = J d'^Xn{X)^{X) - L. Using f(X), Eq. ®, 
the Lagrangian Eq. and using Eq. to eliminate 
V in favor of the conjugate momentum H, to lowest 
(i.e. quadratic) order in v (21 1 we find 



H = Eo[A\+5Hj^[Ii,v\, 
1 



<5i/A = 



dX 



2p 



H g-i H 



Mo 



daV ■ g • daV 



(10) 

(11) 



Here, Eq[A] is the elastic energy for the uniform strained 
reference state, identical to the classical theory result /o, 
Eq. whereas 5H\ describes collective elastic fluctu- 
ations, and depends explicitly on the uniform shear de- 
formation A through the metric tensor g. 

The partition function of a macroscopically sheared 
rubber is then given by the following phase-space path 
integral: 



r/U r/v 



-13SHa 



(12) 
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where (3 = l/fc^T. We note that to quadratic order in 
V and n (which is our focus here) the partition func- 
tion, Eq. ifT^ . separates into a product of kinetic (Z^) 
and elastic (^X) P^^rts. Furthermore, the incompressibil- 
ity constraint only applies to the v field but not to the 
canonical momentum field 11. As a result, the momen- 
tum contribution leads to an inconsequential strain- 
independent constant. We emphasize, however, that be- 
cause of the nonlinear couplings between v and 11, this 
property does not persist to higher orders, and in our 
formulation (in terms of v) the momentum degrees of 
freedom contribute nontrivially to rubber elasticity j^. 
The elastic part of the partition function is given by 



fluctuations: 



7'" 



J DvYl5{q-v) exp -^^y^g^%-g-i7-g 



(13) 



Because of the incompressibility constraint on v, encoded 
in the (5-functional, the dependence of the free energy on 
the imposed strain A cannot simply be eliminated by a 
change of variables. This contrasts with the aforemen- 
tioned result of James and Guth for a phantom net- 
work, where such a constraint, and concomitantly the 
dependence on the imposed deformation, are absent. 

The structure of for the unstrained case of g = I 
is identical to that of a U{1) gauge field theory in the 
transverse gauge. Representing the ^-functional by its 
functional Fourier representation, Z"^ is easily computed 
via two standard Gaussian integrations: 



n I Tr(PL^-.g-i^ 



1/2 



(14) 



where P'" = qq/q^ = qq is the longitudinal projection 
operator onto q. Ignoring irrelevant A-independent ad- 
ditive constants, we obtain the free-energy correction due 
to elastic fluctuations from scales longer than ^: 

SFa - -fcsTlnZX = ifcsr^ln(Trpi^g~i), 



-fcBT^r'f^<i(ln(TrpLg-i))., 



(15) 



where 



denotes an average over the orientation of 



the d-dimensional unit-vector q, while firf is a numerical 
factor of order of the surface area of a d-dimensional unit 
sphere. 

By combining this result with the classical contribu- 
tion, £^o[A] = /o[A], Eq. Q), we arrive at the central 
result of this letter, i.e., the elastic free-energy density 
of an incompressible rubber, subject to a uniform shear 
deformation A, computed to lowest-order in the thermal 



/(A) 



fJ-o 



Trg+-ri7,r'(ln(TrPV)).. (16) 



As was argued earlier, the scale of fluctuation contribu- 
tion (i.e. the second term), measured by Tild£,~'^, is of 
the same order of magnitude as the classical, mean-field 
contribution set by the shear modulus /zq. 

For a three-dimensional uniaxially deformed system, 
with A given by Eq. we have calculated the average 
in Eq. explicitly, and thus obtain an analytical ex- 
pression for /(A). Ignoring a A-independent constant, we 
find 



/ = ^Aio(A' + |)+Mi 



tanh Vl - A- 



- InA 



(17) 



where /ii = 2Tflj^£^^ . For compression, i.e., A < 
1, the second term should be analytically continued 
in such a way that it remains real and positive, 
namely with ( 1 - A'^) '^/^ tanh~^ Vl - A^^ _^ (;)^-3 _ 
l)-i/2 tan-i VA-3 - 1, for A < 1. 

In Fig. n we compare our prediction for the corre- 
sponding Mooney stress, {df/dX)/{\ — A^^), with that 
extracted from the stress-strain curve of Ref. 3], and 
find excellent agreement. The "mean- field" and long- 
wavelength fluctuations shear moduli, fiQ and /ii, re- 
spectively provide two independent fitting parameters, 
which only vertically translate and scale the curve in the 
Mooney-Rivlin plot, Fig. ^ However, the shape of the 
curve is completely determined by the fluctuation con- 
tribution (i.e., second term) of the RHS of Eq. H17(l . and 
thus has no free parameters. Comparison of Eq. H16|l 
with experimental data on biaxial deformations (analy- 
sis of which we leave for the future) should provide more 
stringent test of our theory. 

Our analysis can be easily extended to a more realistic 
system with a finite bulk moduhis B. This can be done 
by removing the hard constraint, Eq. (|SJ| and Eq. (Tj), and 
adding a term B(V • ■i/)^/2 to the Hamiltonian, Eq. Hll|l . 
which suppresses density fluctuations. In the case of a 
uniaxial distortion of a compressible rubber, we find the 
lowest-order fluctuation correction to the free energy be- 
comes 



Sf{X,T) 



Ml 
2 



1- A- 



^/A , -1 
' tanh ^ 



A- 



1 + r/A 



r) 



(18) 



where t = Ho/B. As illustrated in Fig.^ this result nat- 
urally interpolates between an incompressible and phan- 
tom (classical) rubber. This result is consistent with ex- 
periments that observe a systematic reduction in the de- 
viation from the classical theory with softening of the 
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bulk modulus via swelling, i.e. adding solvent to the sys- 
tem [12] . 

The correlation function 

V Gf = {u^ut^) = A„ Ab, {v\vl.) (19) 

of the phonon field u[X) sa A • v(X) (defined by Eq. (jHJ), 
relative to a macroscopically strained state tq = A • AT, 
can also easily be computed. For an incompressible rub- 
ber, a Gaussian integration with the Boltzmann weight 
from Eq. H14|l gives 

^ Mor \ TrP^g ^ J 

There is a small caveat, however. We have been la- 
beling the mass points in the deformed state tq = A ■ X 
by their undeformed equilibrium coordinate X, which is 
conjugate to the wave-vector q. However, the wavevector 
p probed by scattering experiments is the one conjugate 
to the deformed equilibrium position = A ■ X. These 
two vectors are related via 

(p, A • X) = (q, X)^q^A^-p, (21) 

where (•,•) denotes the inner product of two d- 
dimensional vectors. Therefore, the phonon correlation 
function, as a function of the physical wave-vector p is 
given by Eq. (|2()|l with q expressed in term of p through 
Eq. (ini- This leads to 

Gf - ^-..^^^T.p> (^-^ " P'^M « PlbiP)^ (22) 

which is proportional to the transverse projector, as ex- 
pected from incompressibility. The predicted anisotropic 
strain-dependence provides an independent test of the 
theory. 

In this letter we have demonstrated the importance of 
thermal fluctuations for the elasticity of isotropic rubber, 
particularly in the large-deformation regime. It is not dif- 
ficult to see, however, that the same general mechanism 
extends to all incompressible soft solids, such as liquid 
crystalline elastomers. Our current analysis is limited 
to the lowest order. Because the effective coupling con- 
stant is of order unity, we expect high order corrections to 
be quantitatively important. However, as illustrated in 
Fig^ the lowest-order contributions already capture the 
essential effects of thermal fluctuations on rubber elas- 
ticity. Our formalism provides a systematic approach for 
addressing these high order contributions. We hope that 
our work will stimulate further studies in this direction. 
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